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Abstract. Comparison is made between the distribution of saddle points in 
the chaotic analytic function and in the characteristic polynomials of the Ginibre 
ensemble. Realising the logarithmic derivative of these infinite polynomials as the 
electric field of a distribution of coulombic charges at the zeros, a simple mean- 
field electrostatic argument shows that the density of saddles minus zeros falls off 
as 7r _1 |z| — 4 from the origin. This behaviour is expected to be general for finite 
or infinite polynomials with zeros uniformly randomly distributed in the complex 
plane, and which repel quadratically. 



It is well-known that there are several similarities between the distributions of 
zeros of the ensemble of random polynomials which tend to the chaotic analytic 
function (discussed by Hannay 1998, 1996, Bleher, Shiffman and Zelditch 2000, 
Forrester and Honner 1999, Leboeuf 1999, Bogomolny, Bohigas and Leboeuf 1996) 
and the eigenvalues of Ginibre matrices (Ginibre 1965, Mehta 1991 chapter 15), both 
in the finite and infinite cases. To be more precise, for N large and possibly oo, we 
compare the zeros of the chaotic analytic function polynomials (caf polynomials) 



where the a n are independent identically distributed complex circular gaussian random 
variables, with the eigenvalues of matrices in the Ginibre ensemble, defined to be 
N x N matrices with entries independent identically distributed complex circular 
gaussian random variables. The Ginibre analogue to equation (|l|) is the characteristic 
polynomial /jv,Gin(^)- The zeros of the two /jv share the following properties, as 
discussed in the above references: 

• they are uniformly randomly distributed, with density a = 1/ir, within a disk 
centred on the origin of the complex plane, which has radius y/N, and smoothed 
boundary (Ginibre with gaussian tail outside, caf with power law); 

• within this disk, the distribution of zeros is statistically invariant to translation 
and rotation; 

• the statistical properties of the zeros at a fixed radius r do not change as N 
increases, provided that \/N > r + 0(1); 

• the two-point correlation functions g C af,gGin for zeros separated by distance 
\zi — Zi\ = r < \fN —0(1) within the disk are given by 




(1) 



ffcaf (V2r) — [(sinh 2 r 2 + r 4 ) coshr 2 — 2r 2 sinhr 2 ]/ sinh 3 r 2 
ffGin(r) =l-exp(-r 2 ). 



(2) 
(3) 



Letter to the Editor 2 




Figure 1. A sample distribution of zeros and saddles for a caf polynomial with 
N = 100. The zeros are represented by filled circles, the saddles by empty circles. 
The large disk has radius 10, indicating the area in which the distribution is 
isotropic. Lines of constant argument (equally spaced by tt/2) are given by the 
curves; in the electrostatic analogy, they are electric field lines. The plot to the 
right shows the theoretical density of zeros a (full line) and saddles p (dashed 
line) for N = 100, showing the saddle bump at the origin and the disk boundary 
smoothing (the saddle boundary lies just inside the zero boundary, giving one 
fewer saddles than zeros in all. 



Both of these functions exhibit quadratic repulsion at the origin, are of order 1 
(uncorrelated) for r > 1, and satisfy a screening relation: the integral of o~g over 
the plane, after subtracting the uniform background density, is —1; 

• the distribution of Ginibre zeros is equivalent to a two dimensional iV-charge 
Coulomb gas in a harmonic oscillator potential |z| 2 /2, at a temperature 
corresponding to j3 = 2, whereas the caf zeros have an additional TV-body 
potential. 

Despite these similarities, there is no obvious explicit relation between these two 
ensembles. 

In this work, we consider the distribution of the saddle points of the polynomials 
/caf> /Gin (dropping the suffix N unless necessary), that is, zeros of the derivative 
d//dz. The behaviour of the caf saddles is easy to determine, the Ginibre saddles less 
so. Numerical experiment shows that in both cases the saddle distribution roughly 
mimics the zero distribution, except for a surplus near the origin and a deficit near 
the disk edge. Using the electrostatic analogy, we shall see that the density of the 
saddles minus zeros has the same l/7r|,z| 4 tail away from the origin for each case, due 
to the quadratic repulsion of the zeros. Also, the cause and extent of the edge deficit 
will be clear. 
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Figure 2. Plots of Af(r) compared against the asymptotic theoretical fit: (a) caf, 
fitted with 1 — 1/r 2 as in (El) (solid line) and the exact Af from (W), (fc|) (dashed 
line); (b) Ginibre, fitted with 1 — 1/r 2 as in (a). The data in each case are for 
1000 random samples with Af = 50 (numerical errors tend to appear for Ginibre 
with N significantly larger). 



The density of caf saddles p can be found by replacing / ca f with f' cai in the formula 
for the density of zeros (Hannay 1996, Nonnenmacher and Voros 1998 equation (69)). 
When \z\ < \/N — 0(1), the N — > oo formulae may be applied. In this case, the zero 
density a is 1/tt (as mentioned above), and the saddle density is 

p = n~ l d z d z , (fUWUzT) = ^(1 + (1 + \z\ 2 r 2 ). (4) 

This has, in addition to a uniform 1/tt density, a 'bump' in the vicinity of the origin, 
which integrates to 1, and the saddle surplus is p — a = (tt(1 + |z| 2 ) 2 ) -1 , which decays 
like 7r -1 |z|~ 4 . For finite N, of course, an order N polynomial has N — 1 saddles, not 
N + 1, implying that p w for \z\ 2 > N — 2, which is indeed the case when the 
calculation in equation (||) is performed for finite N. The pattern of zeros and saddles 
for a sample caf polynomial with N = 100 is shown in figure [l] (the pattern of Ginibre 
zeros and saddles look very similar), along with the corresponding plots of a and p; 
the most obvious feature of this distribution is that zeros and saddles tend to occur 
in 'bonded' pairs, with the saddle positioned radially inwards from the zero, and the 
bond length decreases as \z\ increases. The pairing breaks down near the origin, indeed 
the 'extra' unbonded zero is in this neighbourhood. 

It proves to be more convenient to work with the cumulative number of saddles 
minus zeros within a radius r = \z\, denoted by Af{r), which is given by 

Af(r) = 2?r / [p - a)r' dr' (5) 
Jo 

w 1 - 1/r 2 for caf, with l«r< \ZjV- 0(1). (6) 

It is our aim to show, using mean-field electrostatics, that this expression holds 
generally for polynomials with zeros distributed according to the bullet points above, 
including the Ginibre polynomials /bin- Figure || shows numerical plots of JV(r) for 
simulations of both the caf and Ginibre polynomials, the caf fitting @ as well as the 
exact theoretical curve (putting the finite form of (||) into (Eh), the Ginibre, fitting 
with (§). 

We start by recalling that the (two dimensional) electric field is represented by 
the logarithmic derivative of the polynomial /, with N zeros Zj , which electrostatically 
represent N unit charges at positions Zj in the plane. The logarithmic derivative £ of 
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/is 

dz z — Zj 

3 = 1 J 

which is Coulomb's law in two dimensions (setting q/2irso = 1), and holds when there 
are arbitrarily many zeros/charges. The direction of the electric vector E is given by 
£* , i.e. E = (Re£, — Im£), and the direction of the field lines is given by the contours 
of constant argument of /. The magnitude of the field \£\ = \E\, and saddles are the 
places where this is zero. It will also be useful to define the field due to all charges 
excluding a particular one at Zi, 

N 1 

*(*)= E —■ ( 8 ) 

3^h3 = l 

Saddles may therefore be thought of as the places where the field from a 
zero/charge balances the background field from all of the others. Consider a particular 
zero at z t with >/~N — O(l) > \z{\ 3> 1. By Gauss's law, the statistically disklike 
distribution of the zeros implies that the average field at zi due to the other charges 
is that due to a charge air\zi\ 2 at the origin (only zeros with modulus less than \zi\ 
are relevant). This justifies electrostatically the observation from figure [I] that, for 
sufficiently large Zi, zeros and saddles are almost always paired, and that the saddle s$ 
paired with the zero Zi is near Zi, on the straight line between the zero and the origin. 
We denote the real positive 'bond length' \z% — S{\ by 6(sj) = bi. This bond length 
function will be calculated below on the basis of the electrostatic model, ignoring 
statistical fluctuations (implying that &(sj) is radial). This will yield Af(r) by finding 
the number #(r) of bonds crossed by a circle of radius r. Since for a polynomial, 
jV(oo) = —1 and #(oo) = 0, we can set Af(r) = #(r) — 1. 

The number #(r) is simply the number of zeros in the annulus whose inner radius 
r is | Si |, and whose outer radius is |sj| + b(si); any zero Zj in this area will have a 
saddle with \sj\ < |sj|, and its bond bj will cross the circle of radius \si\. Therefore, 
with | Si | = r, and realising that the zero density a is uniform, 

Af(r) + 1 = #(r) = ax[(b(r) + rf - r 2 ] 

= b{r) 2 + 2rb(r). (9) 

The problem remains of how to calculate b(r). The crudest approximation is to 
balance radially the field from z\ with the field from the rest of the charges £i, again 
ignoring fluctuations. Using Coulomb's law, and applying Gauss's law at \zi\, gives 
1/6, = \Zi\ or bi = l/\zi\. Solving b t in terms of |sj| and putting into equation ^ does 
indeed give the required leading order terms 1 — 1/r 2 in (JsJ) , but we have made two 
approximations in this argument which affect the value of bi at the required order. 
We show below that these two effects cancel each other. 

The first approximation made is that Gauss's law should really be applied at the 
saddle radius \si\, not the zero radius (\si\ + 6(s,)). Substituting this value into (^|), 
however, gives Af(r) = 1 + 1/r 2 , not 1 — 1/r 2 as desired. The second approximation is 
that the repulsion from Zi, embodied by the correlation function <?, has been neglected. 

The crude approximation assumed that the field due to the charges other than 
Zi to be due to a uniform jellium of density a. However, we know from the two-point 
correlation function g (with the properties above) that zeros are repelled quadratically 
from a given one, and the background jellium is 'dented' around Zf, with the shape 
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of the dent given around Zi by g(\z — Zj\). The correct field to use, in this case, is 
the mean of (||), not ([?]). Gauss's law, now applied to the dent (since it is circularly 
symmetric around Zi), effectively weakens the field from z^. Including this correction 
as well, we have the implicit expression for b (omitting i subscripts): 

7 — ~T~ f (l-9{b'))b'db' = s. (10) 



b b Jo 

Since b is very small (of the order of 1/s), g{b) is proportional to b 2 due to repulsion, 
and the g(b')_p&rt of the integrand may be neglected, integrating to b 2 /2 to leading 



order; thus (10) implies 

\si\+bi \Zi\ 

which is the same as that crudely derived above. Using this mean-field jellium 
approximation, we have therefore justified the numerical observation that Af(r) ~ 
1 — 1/r 2 , (equation (@)) and therefore the density of saddles minus zeros, to leading 
order, decays as l/nr 4 . We make the following observations: 

• The main objection to this derivation, of course, is that statistical fluctuations 
have been neglected throughout the discussion. By Cauchy's theorem, the exact 
saddle excess Af(r) is 27rr(^(r — Zj)~ 2 / J2( r — z o)~ 1 )' w hich we do not know 
how to evaluate. Both numerator and denominator fluctuate violently, and 
approximating the average of the ratio by the ratio of the averages fails; it 
appears that fluctuations deny us information of the density beyond the r~ 4 
leading order term obtained. Incidentally, for less violently fluctuating fields, 
for example from charges on a unit circle instead of a disk (appropriate for 
CUE (Mezzadri, 2002)), this approximation succeeds in reproducing the leading 
behaviour W(r) « 2Trr(J2(r - z^ 2 £(r - z^ 1 *) / {\ £(r - z^ 2 ) = 1/(1 - r) 
for 1/N (1 — r) -C 1 (the ratio has been rationalised since the denominator is 
otherwise zero by Gauss's law). 

• Fluctuations aside, our result follows only by assuming repulsion between the 
zeros, and it is easy to check numerically that for polynomials with zeros 
distributed completely at random in a disk (Poisson distribution), J\f(r) does 
not have the form (|^). The mathematical form of Af(r) is not known for the 
Poisson distribution. 

We remark that in the infinite caf case, the density of saddles minus zeros, using 
equation (0), is uniform on the sphere upon stereographic projection (Needham, 
1997, p 146); the distribution of neither zeros nor saddles is separately uniform. 
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